Superchemistry: dynamics of coupled atomic and molecular Bose-Einstein condensates 
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We analyze the dynamics of a dilute, trapped Bose- 
condensed atomic gas coupled to a diatomic molecular Bose 
gas by coherent Raman transitions. This system is shown to 
result in a new type of 'superchemistry', in which giant col- 
lective oscillations between the atomic and molecular gas can 
occur. The phenomenon is caused by stimulated emission of 
bosonic atoms or molecules into their condensate phases. 
PACS numbers: 03.75.Fi, 05.30. Jp, 03.65. Ge. 

The experimental observation of dilute gas Bose- 
Einstein condensation (BEC) is revolutionizing low- 
temperature physics BEC represents the ultimate 
limit for the cooling of a gas, since nearly all of the atoms 
may occupy the condensate ground state. Rather than 
being identifiable as single particles, these atoms coher- 
ently populate a matter wave field, and can be coupled 
out of the condensate to produce an 'atom laser' This 
suggests the possibility of nonlinear atom optics. Due 
to atomic interactions, even a single-species condensate 
may exhibit nonlinear wave behavior analogous to self- 
phase modulation in optics Q. Here, we propose that 
a more general type of nonlinearity may occur through 
atom-molecule coupling, and show that this coupling may 
result in the formation of a molecular Bose condensate 
through stimulated emission of molecular bosons. 

More generally, we define 'superchemistry' as the co- 
herent stimulation of chemical reactions via macroscopic 
occupation of a quantum state by a bosonic chemical 
species. In other words, 'superchemistry' results in 
greatly enhanced, non-Arrhenius chemical kinetics at 
ultra-low temperatures. In the simplest case oi A + B 
C reactions, there are three possibilities for the quan- 
tum statistics of the components: bb b, bf ^ f, and 
// —^ b, where b stands for bosonic and / for fermionic. 
In all three cases stimulated emission can occur. Inter- 
estingly, the latter two of these cases correspond to well- 
known quantum-field theories, the Lee- Van Hove model 
of meson theory, and the Friedberg-Lee model of high-Tc 
superconductivity . In this Letter we consider a chem- 
ical system of the first type where bosonic enhancement 
of the chemical dynamics is the strongest. We develop 
the theory of coherently interacting atomic and molecu- 
lar condensates needed to describe this process, and con- 
sider a specific coupling mechanism based on stimulated 
free-bound Raman transitions ||] . 

We begin with the usual quantum field theory Hamil- 
tonian for a noninteracting (atomic or molecular) species 
(i), in a well-defined internal state: 
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Here, is the mass, Vi(x) is the trapping potential, and 
Ei the internal energy of species i. The particles also in- 
teract through collisions. We consider particle number- 
conserving collisions mediated by an inter-species poten- 
tial J7y(x), and non-conserving collisions mediated by an 
effective potential Xyfe(x). The first of these nonlinear 
terms describes the well-known intra-species repulsion or 
attraction, as well as inter-species couplings [||. It is 
desirable to introduce a momentum cutoff to simplify 
the field theory |^], and to replace Uij{x) by an effec- 
tive pseudopotential Uij6(x). This describes low-energy 
S'-wave scattering only. Similarly, the potential x can be 
replaced by an equivalent S'-wave pseudopotential, again 
with the proviso that a momentum cutoff is introduced 
at the level of km ^ a^^, where a is the longest scattering 
length in the problem. The result is an effective quantum 
field theory [gj in which: 

^e// = /'^'^ [I>,l(x)$t(x){7,^.<|./X)4,(x) 
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In the diagonal case, Uu = Airh'^ai/mi, where is 
the i-th species scattering length. In the present work, 
we assume that the trap potential Vi(x) is harmonic: 
Vi(x) = (mi/2)cL'^|xp, where uji represents the rota- 
tionally symmetric trap-oscillation frequency for the i-th 
species. 

The new feature introduced here is the particle number 
non-conserving potential x- Terms like this occur in non- 
linear optics, where they describe parametric processes of 
sub- and second-harmonic generation, which change the 
photon number While matter is clearly not created 
or destroyed in low temperature experiments, an anal- 
ogous effect can occur where two atoms are converted 
into one molecule. Inside a Bose condensate, this chemi- 
cal conversion is dominated by coherent stimulated emis- 
sion, in which transitions are enhanced by the number 
of molecules already occupying the ground state. This 
is completely different from the usual chemical kinetics, 
which predicts that the rates of chemical reactions do 
not depend on the number of particles in the product 
mode, and go to zero at low temperatures according to 
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the Arrhenius law. This type of classical (Boltzmann) ki- 
netic theory is inapplicable in BECs, where the particle 
wavelength exceeds the interparticle spacing. 

In general, the conversion process i + j ^ k will be 
non-resonant. The exception, for which Ek — Ei + Ej, 
corresponds to a Feshbach resonance; such resonances 
have recently been studied experimentally [Q. Alter- 
natively, energy conserving transitions are possible if 
Xijk has a harmonic time dependence. In this paper, 
we consider the specific case of stimulated Raman cou- 
pling induced by two laser fields E^i and El2 of fre- 
quencies ujli and ijJl2, as illustrated in Fig. 1. This 
coupling becomes resonant when the Raman detuning 
S = {2Ei ~ E2)/h— (ujl2 — goes to zero. This al- 
lows coupling to a single molecular state, which can be 
selected by the Raman laser frequencies. 

The implication of these new terms is seen most easily 
by considering the corresponding mean field equations, 
in which the operators are replaced by their mean val- 
ues, and a factorization is assumed. Elsewhere using 
a variational technique, we have shown that this gives a 
good estimate of the ground-state energy at high density 
- relative to the spatially uncorrelated behavior that can 
occur at low densities (corresponding to a Bose gas of 
dressed dimers). In the present case the relevant equa- 
tions are obtained, in a rotating frame, for the simplest 
case of one atomic species 0i , together with a correspond- 
ing molecular dimer species (j)2'- 
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Here we assume that U12 = U21 is the only number- 
conserving scattering process, while x = X112 describes 
conversion of atoms into molecules by stimulated Raman 
transitions. Many interesting dynamical properties of 
these types of equations - including nonlinear oscillations, 
non-equilibrium phase transitions, and soliton formation 
- have been explored in nonlinear optics [pl| . A novel 
feature here is the presence of the trap potential which 
localizes the interaction volume. 

We derive the Raman coupling coefficient x for a sim- 
plified model of the two-body interaction, in which the 
atoms interact in their electronic ground state through a 
single Born-Oppenheimer potential Vg{R). Molecules are 
formed in a single bound vibrational state of energy E2 
with radial wave function 7/2 (i?). Two free atoms with 
zero relative kinetic energy have a total energy 2Ei , and 
a relative radial wavefunction ui{R), normalized so that 
asymptotically Ui ^ ^/in{R — ai). We assume that the 
two laser fields E^^ = Eoi cos(ijJLji) = 1,2) couple the 
ground electronic state to a single electronically excited 
state described by a potential Ve (i?) , with Rabi frequen- 
cies = |dM • Eoi\/h where cIm is the electric dipole 
matrix element connecting these two states. The excited 



state has vibrational levels \v) with energy Ey and ra- 
dial wave functions Ue,v{R)- The excited levels decay by 
spontaneous emission at a rate 7m. All bound levels are 
normalized so that J dR\u^^y(2)\^ — 1- 

To proceed further, we first notice that the effective 
Hamiltonian in first-order perturbation theory should re- 
produce the known behavior of two atoms in a relative 
S'-wave scattering process in the presence of an external 
radiation field . Here we recall that the effective field 
theory has a momentum cutoff (otherwise the perturba- 
tion theory would diverge for higher order terms). From 
these requirements we obtain 
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Here Uq/H — Anhai/mi, A„ = {E^ — 2Ei)/h — ojli, 
Ai^^ — {Ey — 2Ei)/h~ujL2, and Ij^y are the overlap inte- 
grals Ij_y — J dRue,y{R)u*{R). In addition, the molecu- 
lar spontaneous emission leads to the incoherent produc- 
tion of molecules in different states, together with atomic 
excited state decays. Treating these as loss processes, we 
obtain additional terms of form: 



= -q;(/)i + i(3i(f)i - Fil^il (1)1, 

4>2 = -r202 + iP24>2, 

where the induced decay rates are: 



(3) r, = 32^V 
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and i(3j4'j is a light shift term, with 
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Here we have introduced AJ, ' = {Ey — E2)/h — ojli, 
Di — ujQ — ujLi, where loq is the resonance frequency of the 
atomic transition between the dissociation limits of the 
excited and ground potentials. Also, = jd^ • Eoi|//i 
is the atomic Rabi frequency, d^^ is the transition dipole 
moment, and is the atomic excited state population 
decay rate. The Raman detuning at trap center for an 
atomic BEC is S = S + P2 - 2(3i + 2{Uii/h)\(t)i{0, 0)^. 

Rotationally or vibrationally inelastic atom-molecule 
collisions may also give rise to losses. The magnitude of 
these rates is presently unknown, and we neglect them 
here. We note that these rates should decrease rapidly 
with increasing molecular binding energy and go to zero 
in the molecular ground state, so that it should be pos- 
sible to obtain a very low rate by selection of an appro- 
priate level. 
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We have calculated Un, x, Fj, and a for a Vg{R) 
which closely approximates the ®'^Rb2 ground poten- 
tial, and a Ve{R) which closely approximates the ®^Rb2 
0~ symmetry potential that connects to the 5^Si/2 + 
5^P3/2 dissociation limit. Free-bound Raman coupling 
of similar states in ^^Rb2 has previously been explored 
experimentally In this calculation, |2) is the bound 
state of Vg corresponding to the vibrational quantum 
number v = 29, with a binding energy of —160.7 GHz 
with respect to 2Ei, and we take 7^ = 3.7x 10 ''s"^, — 
27A, Oi = 2 X 10^° s-\ ^2 = 6.324 x 10^ s-\ and fl^ = 
r2i/V2. We also choose wli = cjq ^ 15.485 cm^^, where 
(jjQ is the resonance frequency of the 5^Si/2 ^ 5^P3/2 
atomic transition. For these parameters, we find that 
x/h ~ 7.6 X 10"^ m3/2/s, Fi ~ 1.629 x lO^^s ni3/s, F2 ~ 
304.4s-\ /3i ~ 2.108X 10^s-\ /32 ~ 3.344x lO^s"!, and 
a = 134.06 s~^. For typical, realizable BEG densities of 
n 4 X 10^" m~3, we find that x"^^^^ > r2,Fin , as 
required for coherent dynamics. 

Based on these results, we have carried out simulations 
ofEqs. P and (|), using a = 5.4 nm ||l|l, C/n/^i ~ 4.96 x 
10"^'' m^s, and with additional terms given by Eqs. (|^). 
The atom-molecule {U12) and molecule-molecule {U22) 
scattering rates are neglected. Provided they are not 
too strong, the effect of these terms will simply be to 
change the condensate self- energies, and hence to modify 
the optimum Raman detuning. We also choose bJi/2'K — 
^02/2,1^ — 100 Hz, and an initial number of atoms TV — 
5 X 10^. We assume an initial condition of a pure atomic 
BEG, with no molecules present, as given by the steady 
state of the standard Gross-Pitaevskii equation in a trap. 

The results are shown in Fig. 2. We observe giant oscil- 
lations between atomic and molecular condensates, which 
take place on short time scales. The integrated atomic 
and molecular numbers, shown in Fig. 3 (a), do not show 
complete atom-molecule conversion, because the oscilla- 
tion frequency is higher in the center of the trap, due 
to the inhomogeneous atomic density. The total number 
shows a small decay due to spontaneous emission. 

The pronounced oscillation between an atomic and a 
molecular condensate provides clear evidence of a long- 
range coherence effect. In contrast, stimulated Raman 
photoassociation in a thermal cloud of atoms would not 
produce similar collective oscillations, because the phases 
associated with the individual atom/molecule conversion 
processes are random in a thermal cloud. The effect is 
also very different from stimulated transitions between 
two spin states of Bose-condensed atoms. That con- 
version is linear in the atomic amplitudes, and therefore 
the stimulated transition rate is equal to the single parti- 
cle rate and independent of the density ||l^. In contrast, 
stimulated atom-molecule conversion in a condensate is 
nonlinear in the atomic and molecular amplitudes, and 
the conversion rate scales with density as fii cx X'^1%/'^- 

Fig. 3 (b) shows the result of the calculation for a re- 
duced atom number and two times smaller density, but 
with the same effective initial Raman detuning. The 
density-dependence of the superchemistry transitions is 
evident in the comparison of Figs. 3(a) and 3(b), which 
shows that the higher density cloud oscillates faster. This 



effect could be studied by optically imaging the atomic 
cloud at a succession of times. Observation of an essen- 
tially pure atomic condensate during one of the atom 
number maxima would imply that the molecular gas 
which exists at a prior time must be coherent. Impor- 
tant signatures of the coherent stimulation are a reaction 
rate that initially increases with time (see insert. Fig. 3 
(a)), and the density dependence of the nonlinear oscil- 
lation period. 

In summary, number-nonconserving interactions be- 
tween bosonic atomic and molecular condensates at 
ultra-low temperatures is predicted to result in a form 
of 'superchemistry' - in which Bose-enhanced coherent 
quantum dynamics replaces the usual chemical kinetics, 
giving a completely new type of behavior. We find that 
giant collective oscillations can occur between atomic and 
molecular Bose condensates. The phenomenon is the 
matter wave analog to optical frequency doubling and 
parametric down conversion. Interesting quantum be- 
havior niay also result, ranging from squeezed-state gen- 
eration ||l^ to quantum soliton formation ||], or even 
quantum chaos. Quantum effects may change short- 
distance correlations, and need to be included in any 
treatment going beyond the mean field theory. In the 
future, the highly specific nature of these proposed Bose- 
enhanced reactions could open the way for new types 
of quantum-controlled chemical synthesis, or other novel 
and unexpected quantum-phase dependent phenomena. 
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FIG. 1. Schematic diagram of the Raman photoassociation. 
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FIG. 2. Densities |(/>f(x,t)|^ for the atomic (a) and molec- 
ular (b) species as depending on time and radial distance 
r = jxj from the trap center, for 5 = 3.879 x lO'^ s~^ (so 
that 5 + /32 - 2/3i = -2.8 x 10^ s"^). 
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FIG. 3. (a) Occupation numbers Ni = J dx\(f)i{x,t)\'^ of 
the atomic (solid line) and molecular (dashed line) fields, as 
a function of time t, for the parameter values of Fig. 2; (b) 
same as in (a) but for the half the initial atomic density and 
the same initial efi'ective detuning 5. 
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